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Abstract Resonant cyclotron scattering (RCS) in pulsar magnetospheres is considered. 
The photon diffusion equation (Kompaneets equation) for RCS is derived. The photon 
system is modeled three dimensionally. Numerical calculations show that there exist not 
only up scattering but also down scattering of RCS, depending on the parameter space. 
RCS's possible applications to the spectra energy distributions of magnetar candidates 
and radio quiet isolated neutron stars (INSs) are point out. The optical/UV excess of INSs 
may caused by the down scattering of RCS. The calculations for RX J1856. 5-3754 and 
RX J0720.4-3 125 are presented and compared with their observational data. In our model, 
the INSs are proposed to be normal neutron stars, although the quark star hypothesis is 
still possible. The low pulsation amplitude of INSs is a natural consequence in the RCS 
model. 

Key words: radiation mechanism: nonthermal, scattering, stars; neutron, pulsars: gen- 
eral, pulsar: individual: RX JI856.5-3754, pulsar: individual: RX J0720.4-3125 



1 INTRODUCTION 

Three kinds of pulsar-like objects have additionally and greatly boosted up our konwledge about pulsar 
magnetospheres. They are anomalous X-ray pulsars and soft gamma-ray repeaters (magnetars candi- 
dates), radio quiet isolated neutron stars (INSs) (the magnificent seven) , and rotating radio transients 
(RRATs). Figure [T] shows their positions on the P — P diagram. Our conventional picture of pulsar 
magnetospheres is provided by e.g. Goldreich & Julian (1969), Ruderman & Sutherland (1975), and 
Cheng et al. (1986) (for a recent review, see Kaspi et al. 2006), which is mainly about the open field line 
regions (OFLRs). Few people begin to realize that there could be interesting physics in the closed field 
line regions (CFLRs) of pulsar magnetospheres. For magnetars, it is proposed that there is strong and 
twisted magnetic field around the central star (Thompson et al. 2002; Lyutikov & Gavriil 2006). INSs 
are thought to be dead neutron stars, which provide a clear specimen for magnetospheric and cooling 
studies (Kaspi et al. 2006; Tong & Peng 2007; Tong et al. 2008). For RRATs, recent modeling also 
indicates interesting physics in CFLRs (Luo & Melrose 2007). The most direct evidence comes from 
the observations of the double pulsar system PSR J0737-3039A/B, and there could be also signatures of 
interesting physics in CFLRs of normal pulsars (Lyutikov 2008). 

The interesting physics in pulsar CFLRs are mainly related to the plasmas there. Roughly speaking, 
the electron density in magnetar CFLRs is about 4-5 orders higher than the Goldreich-Julian density 
(Rea et al. 2008). In the case of RRATs, Luo & Melrose (2007) have proposed an idea of "pulsar 
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Fig. 1 P — P diagram of pulsars. Diamonds are for INSs (Haberl 2007; Kaplan 
& van Kerkwijk 2009). Stars are for RRATs (McLaughlin et al. 2006; McLauchlin 
et al. 2009). Magnetar (squares) and radio pulsar (dots) data are from ATNF 
(jhttp://www. atnf . csiro.au/research/pulsar/psrcat/ 1 . The dotdashed line is the constant poten- 
tial line 1/ = 1/3 X 10^^ V. 



radiation belt", like the radiation belt of the earth. Noting the similarities between INSs and magne- 
tars/RRATs, we suggest that there could also be plasmas in INS CFLRs with number density greatly 
higher than the Goldreich-Julian density (i.e. the electron blanket, see Wang et al. 1998; Ruderman 
2003). We observed three similarities between INSs and magnetars/RRATs: 

1. Most of them are long period pulsars with spin periods about 10 seconds; 

2. They all show a non-atomic Planckian spectrum; 

3. They all have large or relatively large spin-down rates (indicating possible higher fields). 

With these similarities, we suggest that the physics of these three kinds of objects should be similar, and 
that the very different observational manifestations could be resulted from a different parameter space 
and an evolution history. 

In this paper we consider resonant cyclotron scattering (RCS) process in pulsar magnetospheres. 
In section two, a brief description of the RCS process and some basic formulas are given. Kompaneets 
equation for the RCS process is deducted in section three. Numerical calculations are given in sec- 
tion four The application to INSs is the topic of section five. In the last two sections, discussions and 
conclusions are summarized, respectively. 

2 RESONANT CYCLOTRON SCATTERING 

Near the neutron star surface, RCS of photons is more important than Compton scattering (Ruderman 
2003). It has the following three points. 

1 . Given the magnetic field, the scattering occurs only at a specific frequency and vice versa. Given 
the frequency, the scattering occurs only at a specific location in the pulsar magnetosphere. At the 
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resonant frequency, the cross section is about 8 orders larger than the Thomson cross section for 
typical magnetic field 10^^ G. 

2. The momentum is conserved only along the z-direction (direction of the magnetic field). The field 
may absorb or contribute perpendicular momentum. 

3. The particle distribution is strongly affected by the field. The perpendicular motion is suppressed 
while particles can move freely along field lines. Therefore the electron distribution is ID. 

The soft X-ray spectrum of magnetars my be the result of RCS of surface thermal emission (Rea et 
al. 2008). One has three ways dealing with scattering problems (including resonant scattering). One 
is solving the radiation transfer equation directly (e.g. Lyutikov & Gavriil 2006), the second is doing 
Monte Carlo simulations (Fernandez & Thompson 2007; Nobili et al. 2008), and the third is introducing 
photon diffusion equation (Kompaneets equation) as in the Compton scattering case (e.g. Rybicki & 
Lightman 1979). However, the Kompaneets equation for RCS has not been developed yet. Considering 
it's importance in magnetar soft X-ray emission, we present a Kompaneets equation method for the RCS 
process in this paper. Also improved approximations are employed. We find that it may account for the 
optical/UV excess of INSs. 

Previous solution provided by Lyutikov & Gavriil (2006) has three problems which should be im- 
proved. The three problems are: 

1 . It's a one dimensional treatment. Photons can only propagate forward or backward. All the calcula- 
tions and arguments there are valid only in the ID case. This will cause two additional problems. 

2. The angular dependence of the RCS cross section is smeared out. The rigorous expression is eq. (|3]l 
which we will discuss in the following. 

3. The down scattering of photons is dropped. In the ID case, the phase space volume is proportional 
to oc p. While in the 3D case it is proportional to oc p"^. Noting that photons are bosons, the key 
difference between the 3D and ID case is that there is no Bose-Einstein condensation in the later 
case (Pathria 2003). While the down scattering is Bose-Einstein condensation of photons in the low 
energy state (Liu et al. 2004; Sunyaev & Titarchuk 1980), therefore it can not be handled in the ID 
case. Then it is not surprising that the authors found a net up scattering of transmitted flux. It is the 
approximations they employed that matters. We try to provide a 3D treatment of the photon system 
in this paper. 

Before going to deduction details, some basic formulas should be given at first (You et al. 1997). 
Cyclotron frequency of electrons in a given magnetic field is 



where vb is the local cyclotron frequency, lub is the angular frequency lob = ^tti/b, e is the electron 
charge (absolute value), B{r) is the local magnetic field, r is the distance from the point to the center of 
the star, mp is the electron rest mass, and c is the speed of light. When v = vb, where v is the photon 
frequency, RCS occurs. The differential cross section is 
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(1 + cos^ e){l + cos^ e')(j){v - VB)dVL\ 
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where is the classical electron radius, 9 is the angle between the incoming photon and the local 
magnetic field, 0' denotes the angle of the outgoing photon, dVl' is the solid angle of the outgoing 
photon, and (f){v — vb) is the Lorentz line profile function, which acts like a Dirac delta function 
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Note that F is the natural width 



The Lorentz line profile function has the normalization condition 

f + CXD 



— VB)di' = 1. (6) 
Performing the angular integral gives the total cross section 

c^Rcs = ^7rrec(l + cos^ 0)4){v - vb), (7) 

which depends on frequency. 

In the case of pulsars, dipole magnetic field is always a good approximation. The magnetic field at 
radius r is ^ 

B{r)=Bp{^^ , (8) 

where Bp is the magnetic field at the surface of the neutron star, and R is the neutron star radius. Given 
the photon frequency v, the radius at which RCS occurs is 

TKCS = (^) R, (9) 

where is the cyclotron frequency at the star surface i^Bp — 27 ^Tc (oi^ly photons with frequency 
smaller than VBp will encounter RCS). For photons in the soft X-ray band 1 keV < hiy < 10 keV 
(with h is Planck's constant), we are only considering a specific frequency range i^i < v < V2. The 
scattering occurs in a finite space range r2 < r < ri, where r2 is the scattering radius corresponding to 
frequency V2, and ri corresponding to frequency vi. We assume that there are bulk of electrons filling 
the space between r2 and ri . Beyond ri there may also be bulk of electrons, but is is less related to the 
observations in the frequency range vi to V2- Finally we introduce the optical depth of RCS 



trcs = j Necr^csdr = tq{1 



cos^e), (10) 



where is the electron number density (assuming homogenous), 

Tre^TVerRcs 

To = — ^ ■ (11) 

The optical depth also depends on frequency (x l/v^^^. In the following sections, all optical depth are 
referred to their value at the lower frequency boundary, i.e. optical depth at i^i. During the integration of 
eq. ( fTOb . the spatial dependence of magnetic field (eq. ([8])) is taken into consideration. This will be used 
in the numerical calculation section. 



3 KOMPANEETS EQUATION FOR RESONANT CYCLOTRON SCATTERING 

We formulate our deduction analogous to that of Kompaneets equation for Compton scattering (e.g. 
Rybicki & Lightman 1979; You 1998; Padmanabhan 2000). 

Denote the initial and final state of the scattering as {pz, n) and (p^, v\ n'), respectively, with 
Pz is the initial electron momentum in the z-direction, v the initial photon frequency, n the propaga- 
tion direction of the incoming photon, and a prime donotes the corresponding quantity of the outgoing 
particles. In strong magnetic field, electron motions perpendicular to the magnetic field are celled into 
Landau energy levels. Almost all the electrons are in the ground state (You et al. 1997). And we only 
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consider photons in the ground state before and after the scattering (Herold 1979). According to energy- 
momentum conservation in the non-relativistic case, we have 



hv 



2top ' 



(12) 
(13) 



It seems that we are dealing with a ID distribution of electrons. 

From the conservation of energy and momentum, we can calculate the frequency change after and 
before the scattering A = — ly)/ kT^, with k is Boltzmann's constant, and the effect temperature 
of the electron system. Since we are dealing with non-relativistic electrons kT^ WeC^, and consider 
typical photons in the X-ray band hv ^ 1 keV m^c^, the frequency change is very small A <C 1. 
Therefore only considering first order terms of A, we have 



A = (cos 61 -cos 61') 



(14) 



where x is the dimensionless frequency x ~ hv/kTe,- The above expression is accurate to an order 
of 0{:^^^-pjr), which is negligible in the case of magnetars and INSs. The validity of Kompaneets 
equation method in the nonrelativistic case is well established (e.g. eq.(7.53) in Rybicki & Lightman 
1979). 

Let n{iy) denotes the occupation number per photon state of frequency i/. We denote the transition 
probability from an initial state {pz, n) to a final state (p^,, v' ^ n') as dW. Note that the transition 
probability is a microscopic quantity, we always have dW' = dW. Since electrons move freely along the 
z-direction, we describe the electron system as a ID Maxwellian distribution. The number of electrons 
with momentum in the range Pz — Pz + dpz is f{pz)dpz, with 



fipz) = 7Ve(2^mefcre)-l/2e-f^/2™efeT^= 



(15) 



The evolution of the photon spectrum is described by the Boltzmann equation (Rybicki & Lightman 
1979) 



(16) 



where a subscript RCS means the change of occupation number is caused by RCS, n and n' are abbrevi- 
ated forms of n(z^) and n{v'), respectively. For non-relativistic electrons, the frequency change is small 
A ^ 1, we can expand eq. ( fTSI ) to terms of A^ and neglect higher order terms. The change of photon 
occupation number becomes 



dn 

at 



RCS 



dn , 
d-x^-^' 



dpz / dWfipz)^ 



(17) 



-i- W IL Oft . J_ . , 



dpz J dWfipz)A''. 



We denote the two integrals by 



h = 



dpzj dWf{pz)A 
dpzj dWf{pz)A^, 



(18) 
(19) 
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and then eq. ( fTTI i becomes 

dn 

at 



RCS 



dx 



+ n{l + n) 



(20) 



h 



Using the conservation of photon numbers can greatly simplify the subsequent calculations. 

The number of photons is conserved during the scattering process. Thus we have the continuity 
equation of n{x) in the frequency space 



dn 

'at 



(21) 



where j is the photon flux in the frequency space. Assuming n{x) is isotropic (the validity of the 
isotropic assumption will be discussed in the appendix), we have 



dn ^ 9 2 ■\ 



(22) 



A comparison between eq. 
Lightman 1979) 



and eq. ( |20] | gives that the flux j must have the form (Rybicki & 



j{x) = g{x) 



9n , 
-+n(l + n) 



(23) 



Note that in equilibrium conditions, ?i(.t) = (e^ — l)^^- = —n{l + n), we have no "photon flux" in 
the frequency space, j — 0. This is a necessary condition. The same condition can be used to check the 
validity of other forms of diffusion equation (e.g. Liu et al. 2004). Substituting the above equation into 
eq. (|22] | we have 



dn 
'dt 



]_d_ 

x^ dx 



g{x) 



dn 
dx 



+ n(l + n) 



A comparison between the coefficient of in eq. (|24] | and eq. ( |20] | gives g(x) 



(24) 



(25) 



Finally the Kompneets equation for RCS has the form 



dn 
'dt 



J_d_ 

x"^ dx 



x^-h 
2 



dn 
dx 



+ n{l + n) 



(26) 



We only need to calculate the integral I2 ■ 

Substituting the equation of frequency change eq. (fl4] i into the definition of the integral I2 eq. (fT9] l, 
first performing the integral of momentum p^, we obtain 



h = / dWx^N, 



fcTe 



(cos6'-cos6'')^ 



(27) 



The transition probability is directly related to the cross section 

dW = cdo-Rcs 

= c^(l + cos^ e){l + cos^ e')c^{v - VB)dn'. 



(28) 
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Since we are dealing with non-relativistic electrons, the cross section can be approximated by its value 
in the electron rest frame. Performing the integral over dfl', we have 



h = 2x iVeO-RCS c- 



(29) 



where gg is an angle dependent factor ge = h + h cos 9^. Finally the Kompaneets equation for RCS is 



dn 



kTe 1 d 



RCS 



x'^NeCT-Rcs cge 



dn 
dx 



77.(1 + n) 



(30) 



The cross section ctrcs now depends on frequency, therefore it can not be taken out of the curly brackets. 
Except for this difference and an angle dependent factor gg, it is the same as the Kompaneets equation 
for Compton scattering. 

4 NUMERICAL CALCULATIONS 

Before we make numerical calculations of eq. (|30] |. two integrations should be performed. Performing a 
space integral on both sides of eq. ( [3Ql ), employing the concept of optical depth as in eq. ( fTOl i. we obtain 
the Kompaneets equation in the case of pulsars 



dn 

at 



RCS 



nieC^ x^ dx 



4 ■''RCS 

X egg 

n ~ r2 



dn 



77) 



(31) 



The exprsssion ri — r2 in the denominator is the range of integration. It is also the range of cyclotron 
scattering corresponding to the frequency range i/^ to 1^2- The space integration must be employed in 
order to eliminate the Dirac delta function in the RCS cross section. 

There is an angular factor in the RCS optical depth trcS' we denote it as = 1 + cos^ 6. To 
simplify the calculations, we use the average value of fg and gg 



fe = 



98 



(32) 
(33) 



Here 'gg ~ fggg / fg. From eq. dSOl l till now, we have performed two integrations. One is integration over 
space range, the other is averaging over the incoming angle. These two integrations are introduced in 
order to simplify the numerical calculations. 

The Kompaneets equation for RCS is a pure initial value nonlinear partial differential equation. 
Solving the pure initial value problem follows the same routine as the mixed initial value and boundary 
value problem. In the case of Kompaneets equation for RCS, we have to be careful since we are working 
a semi-infinite domain, that is < < 00. In the real case, we are only interested in a finite frequency 
range (e.g. in the cases of magnetars, INSs). Therefore boundary conditions are needed. 

The Kompaneets equation for RCS describes the diffusion of photons in the frequency space. It is 
related to the specific intensity as 

hit) ^ ^^n{iy,t). (34) 



For a blackbody spectrum, the initial condition i£] 



n{x, t = 0) = 



cxp( 



(35) 



' We ai'e considering a spherical shell of electrons, extending from r2 to ri. When the radiation reaches the lower boundary 
r2, this is taken as t = 0. At r2, it is already several radii far away from the neutron star surface. The gravity there is already very 
weak. Therefore general relativistic effect is negligible when the radiation propagates from r2 to r'l . 
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Since we only consider pure scattering between electrons and photons, the number of photons is con- 
served. The choice of boundary conditions must guarantee this requirement. One guess is that there are 
no photons going in or out of the specified frequency range. Mathematically this is 

dnix ti 

— :r-^ \- n{x,t)(l + n{x,t)) ^ 0, when a; = a; i, 2:2 for all t. (36) 

ox 

Compare eq. ( l23T l (Ross et al. 1978). During the numerical calculations, a simplified version is usec0 

dn{x, t) 



dx 



0, when a; = xi, a;2 for all t. (37) 



In obtaining the final RCS modified spectrum, we use the random walk approximation. A photon entered 
the lower boundary r2 escapes the outer boundary ri after an average diffusion time scale 

idif = mcs,[vi){ri - f2)/c. (38) 

In the random walk approximation, n(x, tdit) is the final output. 

Figure |2] and figure [3] are the numerical results for up scattering and down scattering case, respec- 
tively. Figure |2] shows the result for typical parameters of magnetars. During the calculations, the mag- 
netic field and stellar radius are chosen as typical values of B = 4.4 x 10^'* G and R ~ 10® cm. It 
can reproduce a stiffened blackbody spectrum, therefore it may be applied to interpret magnetar soft 
X-ray spectrum. Figure [3] shows the calculation for typical parameters of INSs with B ~ 10^^ G and 
R = 10® cm. It produces a spectrum with optical/UV excess. Therefore it may account for the opti- 
cal/UV excess of INSs. 

The luminosity of photon^ (number of photons per unit time passing a fixed surface) is proportional 
to j^^ x'^n{x, t) dx. We can calculate this integral before and after the scattering to check whether 
the number of photons is conserved during the calculations. In the down scattering case, the specified 
frequency range spans about three orders of magnitude. The number of photons changes less than five 
percent before and after the scattering. While in the up scattering case, the specified frequency range 
only spans one order of magnitude. If we insert boundary conditions at xi and X2, the number of photon 
is not conserved. It is because the specified frequency range is not wide enough. Therefore we insert 
boundary conditions "far away" from the specified frequency range, at a;i/10 and 15 X2- The number of 
photons changes less than one thousandth. 



5 APPLICATION TO ISOLATED NEUTRON STARS 

ROSAT discovered seven radio quiet INSs (Kaspi et al. 2006; Triimper 2005; for recent reviews see 
Haberl 2007; van Kerkwijk & Kaplan 2007). They all show featureless blackbody spectra, with low 
pulsation amplitude, and high X-ray to optical flux ratio. Their spectral energy distributions show that 
many of them have an optical/UV excess with a factor of several (Burwitz et al. 2001; Burwitz et al. 
2003; Motch et al. 2003; Ho et al. 2007; van Kerkwijk & Kaplan 2007). Table[I]shows double blackbody 
fit of RX J1856.5-3754 (J1856 for short) and RX J0720.4-3125 (J0720 for short). In interpreting their 
optical/UV excess, the emission radius is either too small or too big for reasonable star radius. Several 
theoretical models have been proposed (Motch et al. 2003; Ho et al. 2007; Triimper 2005). Considering 
the discrepancy between current theory and observations, we try to provide an alternative one, in which 
the optical/UV excess of INSs may due to magnetospheric processes, i.e. due to down scattering of RCS 
when the surface emission passing through the pulsar magnetosphere. 

^ Take the equilibrium condition for example. n{v) = ^^/^t j — J"- upper frequency boundary hu2 S> fcT^ad. 

n <^ I, the third type boundary condition is reduced to the second type. At the lower boundary hui < kT,.^^, this is a poor 
approximation. However, the number of low energy photons is proportional to x^n{x)dx oc <C 1, which is also of a small 
amount. Therefore we employ the simplified version of boundary conditions. The validity of this simplified approximation is 
discussed at the end of this section. 

^ luminosity=flux X area=7r^ f — ) 47rr^, then integrate over the specified frequency range. 
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Fig. 2 Modified blackbody spectrum due to resonant cyclotron scattering, up scattering case, 
for typical parameters of magnetars. The initial blackbody temperature is 0.2 keV, the electron 
temperature is 10 keV. The solid line is the initial blackbody spectrum. The dashed line is the 
RCS modified spectrum with optical depth trcsI'^i) = 2. The specified frequency range 
is {vi,V2) = (1 keV, 10 keV). The model parameters are {xi,X2) = (0.1,1), (?'2,ri) = 
(8.0,17) X IQScm. 




2eV 5eV lOeV 100eV200eV SOOeV IkeV 2keV 

hv 



Fig. 3 Modified blackbody spectrum due to resonant cyclotron scattering, down scattering 
case, for typical parameters of isolated neutron stars. The initial blackbody temperature is 
70 eV, the electron temperature is 26 eV. The solid line is the initial blackbody spectrum. 
The dashed line is the RCS modified spectrum with optical depth trcs('^i) = 1000. The 
specified frequency range is (z^i, 1/2) = (1.5 eV, 2 keV). The model parameters are {xi, X2) = 
(0.058,77), (r2,ri) = (3.9,42) x 10^ cm. 



FigurelHand figure|5]show the RCS modified blackbody spectrum in the case of J1856 and J0720, 
respectively. It can account for the optical/UV excess in J 1856 and J0720 quite well. We assume that the 
initial spectrum is blackbody. When passing through the pulsar magnetosphere, it is modified by the RCS 
process. Therefore the final observed spectrum is a modified blackbody, with optical/UV excess due to 
down scattering of RCS. Four parameters are needed; the initial body temperature, the temperature of 
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Table 1 Double Blackbody Fit to INS Spectral Energy Distributions. Tx is the high tempera- 
ture component (X-ray), Tq is the low temperature component (optical/UV), seen at infinity. 
i?x and Ro are the corresponding emission radius, seen at infinity. Here all numbers are 
only estimated values, no error bars are given, which are taken from van Kerkwijk & Kaplan 
(2007). 
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eV 
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km 
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J1856 


63 


5.9 


26 


24.3 


161 


J0720 


85.7 


5.7 


35.4 


23.5 


330 



Table 2 RCS Fit to INS Spectral Energy Distributions. Tiad is the initial blackbody temper- 
ature. To is temperature of the electron system, seen at infinity. is the electron number 
density, "norm" is the solid angle of the source seen by the observer Equivalent hydrogen 
column density 71h is the result of wabs model. 









^c 


norm 






cV 


eV 


10^^ cm"-* 




10^" cm-2 


J1856 


61 


26 


1.8 


7r(10km/161pc)^ 


1.6 


J0720 


80 


35.4 


1.8 


7r(10km/330pc)^ 


2.8 



the electron system, the electron number density (assuming homogenous), and a normalization constant. 
We want to point out that our model has the same number of free parameters as the double blackbody fit 
(two temperatures and two normalization constants). The RCS model parameters are given in table|2] 

During the fitting process, the magnetic field and stellar radius are chosen to be typical values 10^'^ G 
and 10 km, respectively (e.g. Haberl 2007). In the case of J1856, the initial blackbody temperature is 
chosen slightly lower than the high temperature component of the double blackbody fit. The temperature 
of the electron system is chosen as the low temperature component of the double blackbody fit and 
kept fixed during the fitting process. The electron number density are assumed to be homogenous. The 
corresponding optical depth is about 1000. The normalization constant is the solid angle of the source 
seen by the observer The photoelectric absorption cross section is from Morrison & McCammon (1983) 
(the wabs model). The hydrogen column density is consistent with previous studies (Burwitz et al. 2003; 
Ho et al. 2007). 

The case of J0720 is similar. This may in part reflect the similarities between these two INSs. The 
optical/UV excess of INSs in our model is due to magnetospheric processes. This means that, in our 
model, the central star can be a normal neutron star, although other possibilities, e.g. a quark star, can't 
be ruled out (Xu 2002, 2003; for a review, see Xu 2009). 

In our model about the INSs, the blackbody spectrum is from the whole stellar surface. When 
passing through the CFLRs of the pulsar magnetosphere, the photons are down scattered by the RCS 
process. It will bring a decrease of high energy photons. This can account for the observed optical/UV 
excess. At the same time, since the X-ray emission is from the whole stellar surface, it can also explain 
the low pulsation amplitude naturally. 

6 DISCUSSIONS 

We consider the RCS process in this paper. Previous work of Ruderman (2003) is mainly qualitative 
consideration. While our paper is a quantitative one and detailed comparison with observational data 
is also presented. Lyutilov & Gavriil (2006) consider the RCS process by solving the radiation transfer 
equation directly. We want to point the differences between our paper and their's. 

1. An improved approximation is used. The photon system is modeled three dimensionally in our 
paper. 
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14 15 16 17 18 

log v/Hz 

Fig. 4 Spectral energy distributions of J 1856. The squares are observational points (only cen- 
tral values are included). The solid line is single blackbody fit to the X-ray data, the parameters 
are given in table [T] (It is not the initial blackbody spectrum in our model, the parameters are 
listed in tablelH) The dashed line is the RCS modified blackbody spectrum, the model param- 
eters are given in table |2] The specified frequency range and model parameters are the same 
as that in figure|3] All observational data are from van kerkwijk & Kaplan (2007). 

2. A different method is employed. While Lyutikov & Gavriil solved the radiation transfer equa- 
tion directly, we employed the Kompaneets equation method. They are independent methods. The 
Kompaneets equation method is much simpler both analytically (compared to solving radiation 
transfer equation directly, this is why we can employ better approximations) and numerically (com- 
pared to doing Monte Carlo simulations). 

3. Different applications are considered. Previous researchers on RCS mainly focus on its application 
to magnetars. We point out that it may also play an important role in the radio quiet INS case. 
Therefore, a magnetospheric model is presented for the optical/UV excess of INSs. 

4. Our calculations show that there exist not only up scattering but also down scattering of RCS pro- 
cess. 

Further Monte Carlo simulations may tell us more about the down scattering of RCS if there exists 
real down scattering of RCS. The approximation of a ID treatment (Lyutikov & Gavriil 2006) may result 
in a negative point of down scattering of RCS. Our deduction is in the 3D case (for the photon system). 
The key difference is the phase space volume (in ID case oc p, in 3D case oc p^). While in eq. ( fTSI l. 
we note a factor (1 + n) appears. This is a pure second quantization effect. The photons are aligned to 
condensate in the low energy sate, and this quantum effect can only play an important role in 3D case. 
Note that there is also no Bose-Einstein condensation in the low dimension case (see Pathria 2003). 

Through this paper we talk about the RCS process in pulsar magnetospheres. We think it may 
be a common process. In different cases, it has different manifestations. In the case of magnetars, we 
have observed a stiffened blackbody spectrum. In the case of INSs, we have observed an optical/UV 
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Fig. 5 Spectral energy distributions of J0720. The optical/UV data of J0720 is more likely 
to be nonthermal. The specified frequency range is the same as that in figure |3] The model 
parameters are {xi,X2) ~ (0.042,56), (r2,ri) = (3.9,42) x 10^ cm. (The space range is 
the same as that in figure [51 since it is only determined by the specified frequency range, see 
eq ®). All observational data are from van kerkwijk & Kaplan (2007). 



excess. These different manifestations can be treated universally using the Kompaneets equation for 
RCS presented in this paper. 

Concerning the magnetospheric properties, now people are thinking about that the CFLRs of pulsar 
magnetospheres are not dead but filled with dense plasma (Ruderman 2003; Luo & Melrose 2007; 
Lyutikov 2008). The plasma can be lO"' — 10^ times denser than the local Goldreich-Julian density. The 
origin of this over dense plasma is the presence of twisted magnetic field lines (in the case of magnetars) 
or magnetic mirroring (in the case of RRATs). When discussing the magnetospheric properties we have 
to be careful. As stated in section two, the scattering radius and the optical depth (or cross section) is 
frequency-dependent. The local Goldreich-Julian density is proportional to the magnetic field. At the 
scattering sphere, from eq.([T]i, it is proportional to the photon frequency. Given that the electron density 
is a constant, the ratio of Ne/uQj varies with frequency as oc l/v.ln the case of magnetars (photon 
energy ranges between 1 keV — 10 kcV), the electron density is 10""* — 10"' times the local Goldreich- 
Julian density. While in the case of INSs, we have a broader frequency range 1 eV — 1 keV and a much 
higher RCS optical depth, about 1000. The corresponding electron density is 10'^ — 10^ times the local 
Goldreich-Julian density. Therefore, a plasma with number density 10^ — 10^ times the local Goldreich- 
Julian density is presented in the CFLRs of magnetars/INSs according to our model. We have computed 
the mass of this dense plasma. Assuming an electron-ion plasma, the total mass is lO'^^ g and lO'^^ g in 
the magnetar case and INS case, respectively. This is consistent with studies in the double pulsar binary 
PSR J0737-3039, Crab giant pulses, and magnetar spectrum modeling (Lyutikov 2008; Rea et al. 2008). 

The presence of dense plasmas in CFLRs of INSs needs further explanations. Unlike the case of 
magnetars, the INSs are believed to be dead NSs (Kaspi et al. 2006; Triimper 2005). For slow rotators 
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like INSs, the magnetic mirroring mechanism comes to work (Luo & Melrose 2007). Therefore, we think 
that the dense plasma in the case of INSs could be due to magnetic mirroring mechanism. The source of 
this dense plasma may be the result of accretion from circumpulsar material, e.g. ISM, fallback disk etc. 
Unlike the case of RRATs, in the case of INSs the radiation belts are not very far away from the neutron 
stars (about 40 stellar radii at the outer edge). We may call it the "inner radiation belt" of a pulsar if we 
call the radiation belt near the light cylinder proposed by Luo & Melrose (2007) the "outer radiation 
belt" of a pulsar Nevertheless, the particle processes are similar The pulsar accretes material from 
the environment which will be accelerated in the "dormant outer gap" (Luo & Melrose). High energy 
curvature photons will collide with surface X-ray photons generating pairs in INS CFLRs. The pair 
plasma will be confined by the magnetic mirroring mechanism. This is the pulsar "inner radiation belt" 
(or "electron blanket", e.g. Ruderman 2003). Semi-quantitative estimates are given in Luo & Melrose 
(2007), Ruderman (2003). It can be as high as IC* — 10^ times the Goldreich- Julian density. The plasma 
is cold since it has undergone a long time of relaxation (INSs are old thermally emitting neutron stars). 
Meanwhile during the scattering process, the photons will push the plasma particles away from the star 
The kinetic energy of particle decreases thus resulting in a low temperature. This also explains why the 
plasma system in INS CFLRs is distributed in a rather wide space range, see caption of figure[3] Similar 
process is also possible in corona of magnetars (Beloborodov & Thompson 2007). 

The last but not least important question is: can a neutron star have a blackbody spectrum which can 
be modified when passing through its magnetosphere? It might not be impossible. The current neutron 
star atmosphere models leave us two questions: one is that a blackbody spectrum fits the observation 
better than that with spectra lines (Ho et al. 2007). The other is that we haven't found a high energy tail 
in INS X-ray spectra (van Kerkwijk & Kaplan 2007). Therefore from the observational point of view, a 
blackbody spectrum is possible. A blackbody-like spectrum could be reproduced in a quark star model 
(Xu 2009). 

7 CONCLUSIONS 

We consider the RCS process in pulsar magnetospheres. The photon diffusion equation (Kompaneets 
equation) for RCS is presented. It can produce not only up scattering but also down scattering depending 
on the parameter space. Its possible applications to magnetar soft X-ray spectrum and INSs are point 
out. 

The application to INSs is calculated in detail. We show that the optical/UV excess of INSs may 
be due to down scattering of RCS. The RCS model has the same number of parameters as the double 
blackbody model. Mean while, it has a clear physical meaning. The initial blackbody spectrum from the 
stellar surface is down scattered by the RCS process when passing through its magnetosphere. This can 
account for the optical/UV excess of INSs. The low pulsation amplitude of INSs is a natural consequence 
in our model. 

The calculations for RX J 1856.5-3754 and RX J0720.4-3125 are presented and compared with their 
observational data. The model parameters for RX J1856. 5-3754 and RX J0720.4-3125 are similar This 
may in part reflect the similarities between these two INSs. Finally, we point out that the quark star 
hypothesis (e.g. Xu 2002) can still not be ruled out. 

The photon diffusion equation (Kompaneets equation) for RCS is calculated semi-analytically. The 
calculations for the magnetar and INS cases are all for surface thermal emission. Of course, its applica- 
tion is not limited to the thermal emission case. 

Using the Kompaneets equation (both the resonant and non-resonant ones, or a unified one which 
will be presented in the future), a thorough and quantitative study of scattering process in pulsar mag- 
netospheres could be possible. This can help us make clear the physical process in CFLRs of pulsar 
magnetospheres. 
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Appendix A: PROPAGATION EFFECT 

The validity of isotropic assumption employed in the main body is discussed. From eq.(l2T]) to eq.(l22]l. 
we employ the isotropic assumption. It is also implicitly assumed during the angular average in the 
numerical calculation section. Its validity is acceptable in regions not far away from the star This is just 
the case at the inner edge r2 of the "electron blanket". However, at the outer edge ri, its validity needs 
further confirmation. Our approach to this problem is that we consider an isotropization process. From 
eq. (O, the angular dependence of the outgoing photons is (1 + cos'^ 9'), the same as that of cyclotron 
radiation. It is almost isotropic. Therefore, weakly dependent of the angle of incoming photons, the 
photons become isotropic through the RCS process. This provides the required isotropic photon field to 
be scattered by electrons nearby. The upper isotropization process is valid from one space location to 
another. Therefore, the isotropic assumption is valid through the whole space range, from r2 to ri . 

Noting that the number of photons are conserved during the scattering process. The r^^ dependence 
of the solid angle of the star is "transformed" to the photon occupation number n{x,t). It will only 
modify the space integrated form of Kompaneets equation, eq.(l3T]). However, since the output is not 
sensitive to where we introduce the dependence, the results should be similar Detailed calculation 
is presented in below. 



The energy density of the radiation field is 



47r 

— J., (A.l) 

c 



Ju^-^ I ludn, (A.2) 



where is the mean intensity 

J.. = 

47r 

is the specific intensity. From the definition of photon occupation number, we have 

Ui, = nlv) — 5— n.!^, (A. 3) 

where n{v) is the photon occupation number in the Kompaneets equation. Combining eq. ( lA.lb and 
( IA.3I ), we obtain the relation between mean intensity and occupation number 

J, = ^n{v). (A.4) 

For an isotropic radiation field (e.g. as we have assumed in the main text), this is just eq.(l34b. We 
consider the propagation effect for a uniformly bright sphere with brightness By and radius R (e.g., 
Rybicki & Lightman 1979, section 1.3). The mean intensity at radius r is 

Ju - \By{l - Vl - (A.5) 



The dilution factor is 



r 



Jy{r = R) 2 
This is also the dilution factor for the occupation number n{v) 



(A.6) 
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Only considering the propagation effect, it will introduce a spatial dependence of the photon occu- 
pation number 

d 

—r^n(v,r)^0. (A.7) 
or 

For neutron star with surface temperature Trad, the photon occupation at radius r is 

1 1 /i?^ ^ 



Besides the dilution factor, it means that only half of the photons will propagate towards the observer 

The photon occupation number now in eq. dSOl l is a function of frequency, time, and position n = 
n{v^ t, r). In order to include the dilution effect, we introduce another variable m 

m(i^, t) = r'^n{i^, t, r). (A.9) 

From eq. (IA.7l l. m only depends on frequency and time. Multiple on both sides of eq. ( |30] | and 
performing a spatial integral from 7'2 to ri, we obtain 



dm \ _ kTe 1 d f 4 trcs 
/ RCS meC^ dx\ r\— 



dm ,^ m 
dx 



+ m(l + ■ 2 



(A. 10) 



It is similar to eq. (ISTT i. 

In order to make a comparison with the observational data, the flux of such a system must be 
calculate. It is related to the energy density 

= u^c (A. 11) 

= ^n(y,t,r) 

STrhv^ m^u, t) 

Here D is distance of this source. The spectrum is proportional to v^m{v, t). Similar results are obtained 
as in the main text, with similar input parameters. 

There is a deep reason why the isotropic assumption is still valid in regions far away from the star 
From eq. (|9]), low energy photons will be scattered in outer regions. At the same time, they have large 
cross section and optical depth, e.g. see eq. ( fTOl i and the caption of figure [3] Low energy photons will 
encounter strong scattering, although they are limited in a narrow beam. Related issue has already been 
pointed out in section two. 
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